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Abstract: We survey recent work on the properties and decays of mesons as bound states of con- 
fined quarks at both zero and finite temperature. The framework for these investigations is the set 
of QCD Dyson-Schwinger equations truncated to ladder-rainbow level. The infrared structure of the 
ladder-rainbow kernel is described by two parameters; the ultraviolet behavior is fixed by the one-loop 
renormalizaton group behavior of QCD. The work is restricted to the u, d and s quark sector and allows a 
Poincare-covariant study of the masses and electroweak decay constants of the ground state pseudoscalars 
and vectors: 7r, K , p/u>, K* and 4>. Within the impulse approximation, we summarize results for the 7r 
and K charge form factors. Their timelike behavior exhibits the vector meson production resonances and 
from this the associated vector meson strong decay constants are extracted. The confining property of 
the quark self-energy dynamically produced in this work prevents the appearance of spurious qq produc- 
tion effects. The finite temperature behavior of this model, and closely related ones, is summarized for: 
the chiral restoration transition, the deconfinement transition, the T-dependence of the spatial meson 
mode masses, and the T-dependence of selected electroweak and strong decay widths. The T-dependence 
obtained for meson masses is compared to results from lattice QCD. 
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1. Introduction 

The light pseudoscalar mesons play an important 
role in understanding low-energy QCD. They are 
the lightest observable hadronic bound states of 
a quark and an antiquark, and are the Goldstone 
bosons associated with chiral symmetry breaking. 
Their static properties such as the mass and decay 
constants have been studied extensively and are 
qualitatively understood within QCD. A number 
of studies of pseudoscalar mesons [1-3] have played 
a key role in the development of continuum meth- 
ods for modeling QCD via the Dyson-Schwinger 
equations [DSEs]. Such methods have now evolved 
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into an excellent tool for the study of nonpertur- 
bative aspects of a variety hadronic properties and 
processes in QCD [4-6]. 

Here we recall recent studies of light mesons, 
both at zero temperature and at finite tempera- 
ture, in which the quark propagator DSE in the 
it, d and s quark sector is coupled to the meson 
bound state Bethe-Salpeter equation [BSE]. The 
bound state amplitudes so produced are used to 
study meson decays and form factors. This ap- 
proach is consistent with quark and gluon confine- 
ment [5-7], generates dynamical chiral symmetry 
breaking [8], and is explicitly Poincare covariant. 
It is straightforward to implement the correct one- 
loop renormalization group behavior of QCD [3], 
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and to obtain agreement with perturbation the- 
ory in the perturbative region. Provided that the 
relevant Ward-Takahashi identities [WTI] are pre- 
served in the truncation of the DSEs, the corre- 
sponding currents are conserved. The axial WTI 
ensures the Goldstone nature of the pions [9] ; the 
vector WTI is important for electromagnetic cur- 
rent conservation. 

We show results for the pion and kaon electro- 
magnetic form factors in impulse approxima- 
tion [10,11], both in the timelike and the space- 
like region. In impulse approximation, with 
quark propagators dressed in rainbow approxima- 
tion and meson Bethe-Salpeter amplitudes pro- 
duced in ladder approximation, the vector WTI 
for the quark-photon vertex ensures that the me- 
son electromagnetic current is conserved. A fea- 
ture of this work is that the confining property of 
the dynamically generated quark self-energy pre- 
vents the appearance of spurious qq production ef- 
fects in the timelike region. The ground state vec- 
tor mesons provide the lowest physical resonances 
in the structure of the timelike charge form fac- 
tors of the pion and kaon. We summarize the 
ladder-rainbow DSE-BSE results for the ground 
state vector mesons [12] as well as their strong 
and electroweak decays. The various strong and 
electromagnetic interactions involving vector and 
pseudoscalar mesons have been a typical testing 
ground for non-perturbative QCD models [4] and 
here we give a few examples of modern treatments 
of such processes. 

At finite temperature we discuss chiral res- 
toration and deconfinement using the Matsubara 
formalism. In rainbow-ladder truncation, there 
is a second-order, mean field chiral phase transi- 
tion, with a critical temperature around 120 ~ 150 
MeV depending on the details of the ladder ker- 
nel [13]. We also summarize results for the tem- 
perature dependence of the spatial masses of the 
scalar, pseudoscalar, and vector qq bound states 
and selected electroweak and strong decay widths 
of these spatial qq modes [14,15]. Despite decon- 
finement, these modes seem to persist above the 
chiral and deconfinement phase transition. The re- 
sults for the T-dependence of spatial meson masses 



are in reasonable qualitative agreement with re- 
sults from lattice QCD. 

In the next Section we give the relevant DSEs, 
followed by a discussion of the truncation and the 
model parameters in Sec. |3[ In Sec. |we discuss in 
detail the application of the DSE approach to cal- 
culate meson electromagnetic form factors. Modi- 
fications of the model to study QCD at finite tem- 
perature are considered in Sec. [5] and there we also 
discuss the chiral phase transition. Sec. ^ covers 
the temperature dependence of properties of spa- 
tial qq modes, including masses and decays. Fi- 
nally we conclude in Sec. [?]. 

2. Dyson— Schwinger equations 

The dressed quark propagator S(q) is the solution 
of the DSE 1 

Sipy 1 = Z 2 ip' + Z4,m(fj,) 
+ Z l / g 2 D^(k)^S(q)Tl(q,p), (2.1) 

Jq 

where D fllJ {k = p — q) is the renormalizcd dressed- 
gluon propagator, and T l u (q,p) is the renormal- 
ized dressed quark-gluon vertex. The notation 
/ = J d 4 q/(2Tr) 4 stands for a translationally in- 
variant regularization of the integral, with A being 
the regularization mass-scale. The regularization 
can be removed at the end of all calculations, and 
after renormalization, by taking the limit A — ► 00. 
The solution of Eq. ( |2.l[ ) is renormalized accord- 
ing to S(p)^ 1 = m{p) at a sufficiently large 
spacelike /1 2 , with m(/x) the renormalized quark 
mass at the scale /i. In Eq. (|2.l|), S, T 1 ^ and m(/x) 
depend on the quark flavor, although we have not 
indicated this explicitly. The renormalization con- 
stants Z2 and Z4 depend on the renormalization 
point and the regularization mass-scale, but not 
on flavor: in our analysis we employ a flavor- 
independent renormalization scheme. 

The Bethe-Salpeter amplitude [BSA] for a 
meson bound state of a quark a having momen- 
tum g + and an antiquark b having momentum 

1 We use Euclidean metric {7^,71/} = 2<5 M „, 7^ = 7^ and 
a ■ b = Y^ 4 , (lib;. 
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is denoted by T ab (q + ,q^). The meson momen- 
tum is P = q+ — q- and satisfies P 2 — — m 2 . The 
BSAs are solutions of the homogeneous BSE 



K(p,q;P) X a \q + ,q-), (2.2) 



which has solutions at discrete values of P 2 only. 
Here, K is the renormalized qq scattering kernel 
that is irreducible with respect to a pair of qq 
lines, and x a \q,q') = S a {q)T ab (q, q')S b {q') is the 
Bethe-Salpeter wavefunction. In order to imple- 



ment Eq. ( p.2[ ) while preserving the constraint P = 
q+ — q-, it is necessary to specify how the total 
momentum is partitioned between the quark and 
the antiquark. The general choice is q+ = q + r\P 
and q_ = q — (1 — r\)P where r\ is the momentum 
partitioning parameter. The choice of n corre- 
sponds to the choice of relative momentum; phys- 
ical observables should not depend on the choice. 
This provides us with a convenient check on the 
accuracy of the numerics in an actual model cal- 
culation. 

The meson BSA T ab is normalized according 
to the canonical normalization condition 



2 P u = N c -J?- 

dP u 



tr B [f ab (q\ q) S a (q + ) T ab (q, q') S b (q^)] 



k.q 



tr s [x ab Ck' 1 ~k)K(k 1 q- 1 P) X ah (q,q')] 



(2.3) 



at the meson mass shell P 2 = Q 2 = — m 2 , with 
q = q + rjQ, q' = q — (1 — r))Q, and similarly for 
k and k! . The notation tr s denotes a trace over 
Dirac spin. For the antimeson BSA we have used 
the notation T(q_,q + ) = [C~ 1 T(-q + , -qJ)Cf in 
which C — 7274 is the charge conjugation matrix, 
and X 1 denotes the matrix transpose of X. 

For pseudoscalar bound states the BSA is 
commonly decomposed into [3] 

T(fc + r)P, k - (1 - rj)P) = 
l5 [iE+ PF+ jtG + a^ k^P v H] , (2.4) 



where the invariant amplitudes E, F, G and H 
are Lorentz scalar functions f(k 2 ,k ■ P',rj). The 
dependence of the amplitudes upon k ■ P can be 
conveniently represented by the following expan- 
sion based on Chebyshev polynomials 

oo 

f(k 2 ,k p)=J2 Ui(cose) (kpy fi(k 2 ) , (2.5) 

i=0 

where cos 8 = k ■ P j (kP). For charge-parity eigen- 
states such as the pion, each amplitude E, F, G, 
and H will have a well-defined parity in the vari- 
able k ■ P if one chooses r\ = |. In this case, these 
amplitudes are either entirely even (E, F, and H) 
or odd (G) in k ■ P, and only the even (E, F, and 
H) or odd (G) Chebyshev moments /, are needed 
for a complete description. 

3. Ladder-rainbow model 

We employ a model that has been developed re- 
cently for an efficient description of the masses and 
decay constants of the light pseudoscalar and vec- 
tor mesons [3,12]. This consists of the rainbow 
truncation of the DSE for the quark propagator 
and the ladder truncation of the BSE for the me- 
son BSAs. The required effective qq interaction is 
constrained by perturbative QCD in the ultravio- 
let and has a phenomenological infrared behavior. 
In particular, the rainbow truncation of the quark 



DSE, Eq. (2.1), is 



Z X g l D, lv {k) Vl(q,p) - G{k z ) Df^(k) 

(3.1) 

where D l ™°(k = p — q) is the free gluon propagator 
in Landau gauge. The consistent ladder truncat- 



ion of the BSE, Eq. fl2.2|), is 

K(p, q- P) -> ~9{k 2 ) D^ c (k)^ ® f lv , (3.2) 

where k = p — q. These two truncations are con- 
sistent in the sense that the combination produces 
vector and axial-vector vertices satisfying the re- 
spective WTI. In the axial case, this ensures that 
in the chiral limit the ground state pseudoscalar 
mesons are the massless Goldstone bosons associ- 
ated with chiral symmetry breaking; with nonzero 
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current quark masses it leads to a generalisation 
of the Gell-Man-Oaks-Renner relation [9]. In the 
vector case, this ensures conservation of the meson 
electromagnetic current in impulse approximation 
as described in the next section. 

The model is completely specified once a form 
is chosen for the "effective coupling" Q(k ). The 
ultraviolet behavior is chosen to be that of the 
QCD running coupling a(k 2 ); the ladder-rainbow 
truncation then generates the correct perturba- 
tive QCD structure of the DSE-BSE system of 
equations. The phenomenological infrared form 
of G(k 2 ) is chosen so that the DSE kernel contains 
sufficient infrared enhancement to produce an em- 
pirically acceptable amount of dynamical chiral 
symmetry breaking as represented by the chiral 
condensate [16]. 

We employ the Ansatz found to be successful 
for the light mesons in Refs. [3,12] 

4tt 2 D k 2 



Q(k 2 ) 
k 2 



4n 2 lm T(k 2 ) 



hi 



k 2 /A 2 



QCD 



, (3.3) 



with 7 m = 



12 



and F( S ) = (1 - cxp(-^))/ S . 



The first term implements the strong infrared en- 
hancement in the region < k 2 < 1 GeV 2 required 
for sufficient dynamical chiral symmetry breaking. 
The second term serves to preserve the one-loop 
renormalization group behavior of QCD. We use 
m t = 0.5 GeV, t = e 2 — 1, Nf = 4, and we take 
Aqcd = 0.234 GeV. The renormalization scale is 
chosen to be p = 19 GeV which is well into the do- 
main where one-loop perturbative behavior is ap- 
propriate [3,12]. The remaining parameters, lo = 
0.4 GeV and D = 0.93 GeV 2 along with the quark 
masses, are fitted to give a good description of 
(qq), nrin/K and /„•. The subsequent values for fx 
and the masses and decay constants of the vector 
mesons p, </>, K* are found to be within 10% of the 
experimental data [12], see Table [l| A detailed 
analysis of the relationship between QCD and this 
Landau gauge, rainbow-ladder truncation of the 
DSEs with renormalization group improvement, 
can be found in the originating work [3]. A com- 





experiment 


calculated 




(estimates) 


(t fitted) 


rn u=d 
''V=lGoV 


5 - 10 MeV 


5.5 MeV 


m t=lGeV 


100 - 300 MeV 


125 MeV 


- (qq)l 


(0.236 GeV) 3 


(0.24lt) 3 




0.1385 GeV 


0.138+ 


U 


0.131 GeV 


0.131+ 


m K 


0.496 GeV 


0.497+ 


Ik 


0.160 GeV 


0.155 


m p 


0.770 GeV 


0.742 


f P 


0.216 GeV 


0.207 




0.892 GeV 


0.936 




0.225 GeV 


0.241 




1.020 GeV 


1.072 


u 


0.236 GeV 


0.259 



Table 1: Overview of the results of the model for 
the meson masses and decay constant, adapted from 
Ref. [12]. The experimental value for the condensate 
is taken from Ref. [20], the other experimental data 
are from Ref. [21]. 

parison with lattice QCD results [17] confirms the 
qualitative behavior for the dressed quark propa- 
gator produced by this model [18,19]. 

Recent reviews [5,6,22] put this model in a 
wider perspective. These reviews include a compi- 
lation of results for both meson and baryon physics 
with similar models, an analysis of how quark con- 
finement is manifest in solutions of the DSEs, and 
both finite temperature and finite density exten- 
sions. The question of the accuracy of the ladder- 
rainbow truncation has also received some atten- 
tion; it was found to be particularly suitable for 
the flavor octet pseudoscalar mesons and for the 
vector mesons, since the next-order contributions 
in a quark-gluon skeleton graph expansion have a 
significant amount of cancellation between repul- 
sive and attractive corrections [23]. 

4. Electromagnetic form factors in im- 
pulse approximation 

The 3-point function describing the coupling of a 
photon with momentum Q to a pseudoscalar me- 
son, with initial and final momenta P ± Q/2, can 
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be written as the sum of two terms 

Af(P, Q) = Q a A aba + Q b A abb , (4.1) 

where Q is the quark or antiquark electric charge, 
and where A aba (P, Q) and A abb (P, Q) describe the 
coupling of a photon to the quark (a) and anti- 
quark (b) respectively. The meson form factor is 
defined by 



A a u b (P,Q)=2P u F(Q'), 



(4.2) 



with the charge radius defined by r 2 = —6F'(0). 
Analogously, we can define a form factor for each 



of the two terms on the RHS of Eq. (|4.l|), for ex- 
ample 

Af\P,Q) = 2P u F aU (Q 2 ). (4.3) 

Current conservation dictates that each of the 
form factors F a ii{Q 2 ) and F a i a {Q 2 ) is one at 
Q 2 = 0. 

The impulse approximation allows form fac- 
tors to be described in terms of dressed quark 
propagators, bound state BSAs, and the dressed 
qq-f- vertex. We denote by T° (q+ , g_ ) the vertex 
describing the coupling of a photon with momen- 
tum Q = q + — g_ to a quark of flavor a with final 
and initial momenta q+ and g„ respectively. With 



this notation, the vertices in Eq. (4.1) take the 
form [11] 

Af\P, Q) = N c f\ [S a (q) T ah (q, q+) S b (q+) 

Jk 

xiTl(q + , q _)S b (q^T a \q^q)] , (4.4) 

where q = k + P/2 and q± = k - P/2 ± Q/2. The 
expression for A^ ba is analogous. 

4.1 The quark-photon vertex 

The quark-photon vertex is the solution of the 
renormalized inhomogeneous BSE with the same 
kernel K as the homogeneous BSE for meson 
bound states. That is, for photon momentum Q, 

r^(p+,p_) = Z 2l ^+ [ K(p,q;Q)S a (q+) 

Jq 

xr;(? + , ? -)5"(4-), (4.5) 



with p± — p ± Q/2 and q± = q± Q/2. Because of 
gauge invariance, it satisfies the WTI 



,r a j P+ ,p_) = s-\p + )-s- 1 (p- 



(4.6) 



The general form of the quark-photon vertex 
r^(q;Q) requires a decomposition into twelve in- 
dependent Lorentz covariants, which can be made 
from the three vectors 7^, the relative momentum 
q^, and the photon momentum Q M , each multi- 
plied by one of the four independent matrices 1, ^, 
Q, and a^q^Qjj. Four of the covariants represent 
the longitudinal components which are completely 
specified by the WTI in terms of the quark prop- 
agator and they do not contribute to elastic form 
factors. The solution of the BSE for the trans- 
verse vertex can be expanded in eight covariants 
T^(q; Q) with the corresponding amplitudes being 
Lorentz scalar functions. The choice of covariant 
basis T^(q; Q) is constrained by the required prop- 
erties under Lorentz and CPT transformations, 
but is not unique. A convenient basis that facil- 
itates projection of the BSE, Eq. ( |4.5| ), has been 
given in Ref. [10] and a discussion of the results 
for the vertex can also be found there. 

Note that solutions of the homogeneous ver- 



sion of Eq. flh5|) define vector meson bound states 



with masses my = — Q 2 at discrete timelike mo- 
menta Q 2 . It follows that has poles at those 
locations and, in the vicinity of the bound states, 
behaves like [10] 



Yf v {p + ,p.)f v m v 



Q 2 +m 2 / 



(4.7) 



where T^ a v is the aa vector meson BSA, and fy 
is the electroweak decay constant. For the photon 
coupled to u- and d-quarks, this results in a p- 
meson pole at Q 2 = — 0.6GeV 2 . We will come 
back to this behavior in the timelike region in 



Sec. 4.5 



4.2 Current conservation 

At Q = the quark-photon vertex is completely 
specified by the differential Ward identity 



d 



(4.8) 
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If this is inserted in Eq. (4.4), one finds after a 
change of integration variables k — » q — P/2 



A? h (P,0)=2P ll F aU (0) 



N c I tr« 
'q 



T a \q',q)S a (q)T ab (q,q') 



f\ dS»(q>) 



dP„ 



(4.9) 



with q' = q — P. Comparing this expression to 
Eq. (|2]) with 7] — 0, we recognize that the physi- 
cal result F(Q 2 = 0) = 4 follows directly from the 
canonical normalization condition for r afc if the 
BSE kernel K is independent of the meson mo- 
mentum P and the quark-photon vertex satisfies 
the differential Ward identity [24] . 

With an arbitrary value for the momentum 
partitioning parameter rj, the relation between the 
normalization condition for r afc and electromagn- 
etic current conservation is not so straightforward. 
A shift in the value of rj in loop diagrams (with- 
out external quark lines) is equivalent to a shift 
in integration variables. For processes that are 
not anomalous, loop integrals are independent of 
a shift in integration variables provided that all 
approximated quantities in the integrand respect 
Poincare covariance. In this respect, particular at- 
tention must be paid to the representation and ap- 
proximation used for the BSAs. The general BSE 
vertex solution T(q+,q-), considered as a func- 
tion of the incoming and outgoing quark momenta, 
does not depend on 11. It is only in commonly used 



decompositions such as Eq. (2.4), that the value 
of -q becomes relevant. The amplitudes E, F, G, 
and H are scalar functions of k 2 and k ■ P which 
do depend on rj, i.e., on the choice of a relative 
momentum. Under a change of rj, some of the dif- 
ferent Dirac covariants and associated amplitudes 
(such as F and G) will mix, as will the various 
Chebyshev moments /j in Eq. ( |2.5p . However, the 
net result from use of a complete representation of 
both the Dirac matrix structure and the k ■ P vari- 
able is independent of r\ and choice of integration 
variables in loop integrals; this has been checked 
explicitly for the decay constants [3] and for form 
factors [44]. 

The ladder BSE kernel is independent of the 




Figure 1: The two leading-order vertex corrections to 
the rainbow DSE (top) and the corresponding five dia- 
grams to be added to the ladder BSE kernel (bottom) 
for consistency with the relevant WTIs. The quark 
and gluon lines indicate dressed propagators in this 
and the subsequent figures. 



meson momentum P, and with quark propagators 
dressed in rainbow approximation and the dressed 
quark-photon vertex produced in ladder approx- 
imation, the WT4 is satisfied [40]; thus impulse 
approximation for the charge form factor in com- 
bination with rainbow- ladder truncation is ideal in 
the sense that the resulting meson electromagnetic 
current is conserved. That is, the correct electric 
charge of the meson is produced independent of 
the model parameters, and independent of the mo- 
mentum partitioning parameter rj. 

4.3 Beyond rainbow-ladder truncation 

If one goes beyond the rainbow-ladder truncation 
for the DSEs for the propagators, BSAs and quark- 
photon vertex, one has to go beyond impulse ap- 
proximation for the form factors to ensure current 
conservation [44]. For example, one could include 
corrections to the rainbow-ladder DSE and BSE 
kernels that are higher-order in a s , as depicted 
in Fig. 0. Following the general procedure de- 
veloped in Ref. [23], one can show that both the 



WT4, Eq. fl4.6D, and the differential Ward identity, 
Eq. (4.8), are preserved in the truncation indicated 



in Fig. [I]. Also preserved is the axial- vector WT4, 
which is important for the Goldstone nature of the 
pions [9]. 

The resulting BSE kernel K(q,p;P) now be- 
comes dependent on the meson momentum P, 
which means that the second term of the normal- 



ization condition, Eq.(2.3), is nonzero. To be spe- 
cific, with the choice 77 — 0, this introduces the 
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Figure 2: The four diagrams arising from the P- 
dependence of the kernel in the normalization condi- 
tion, Eq. (2.3), if one includes the diagrams of Fig. |l| 
into the DSE dynamics, and chooses all of the meson 
momentum P to flow through one quark only. The 
derivatives with respect to P are marked by slashes. 




Figure 3: The four corrections to the impulse ap- 



proximation, Eq. (1.4), necessary to maintain current 
conservation if one includes the diagrams of Fig. [l] into 
the DSE dynamics. 



four extra terms in the normalization condition de- 
picted in Fig. H These four additional diagrams 
are generated from the BSE kernel in the bottom 
part of Fig. [l] by taking the derivative with respect 
to the meson momentum P, where P flows through 
one quark propagator only. Since a derivative 
with respect to P is equivalent to the insertion of a 
zero-momentum photon according to the differen- 



tial WTI, Eq. (4.6), it is obvious which diagrams 
have to be added to the impulse approximation 
for the vertex A abb to maintain current conserva- 
tion [11], and they are displayed in Fig. ||. In the 
limit Q — > these four additional diagrams be- 
come identical to the four additional diagrams in 
Fig. ||, provided that the vertex satisfies the differ- 
ential WTI. Of course, there are similar contribu- 
tions to A aba , which can be identified with terms 
in the normalization condition with i] = 1. 

A related topic is the correction to the impulse 
approximation from pion and kaon loops. Simple 
addition of such loops without supplementing the 
ladder-rainbow truncation for the DSEs, will gen- 



erally violate current conservation. A consistent 
treatment of the kernels for both the DSE and BSE 
equations and the approximation for the photon- 
hadron coupling is necessary for current conserva- 
tion. At present it is not clear how to incorporate 
meson loops self-consistently in such an approach, 
but we expect corrections coming from such loops 
to be small in the spacelike region. In Ref. [25] it 
was demonstrated that the dressed quark core can 
generate most of the pion charge radius, and that 
pion loops contribute less than 15% to r 2 . For 
larger values of Q 2 the effect from meson loops re- 
duces even further, and for Q 2 > 1 GeV 2 we expect 
the contribution of such loops to be negligible. 

4.4 Results for the tt and i"T ± form factors 

The pion and kaon form factors are given by 



FkAQ 2 ) 
f k °(Q 2 ) 



2 -\F uSu {Q 2 ) + \F uM {Q 2 ), (4.10) 



^F uSu (Q 2 ) + -F uS rs(Q 2 ) , (4.11) 
-\F dM {Q 2 ) + \F d - s - s {Q 2 ) , (4.12) 



where the quark and antiquark charges are evi- 
dent. It is only recently that a consistent imple- 
mentation of the impulse approximation for these 
form factors has been carried out with solutions 
of the ladder-rainbow DSEs for each required ele- 
ment: propagators, meson BSAs, and the photon- 
quark vertex [11,10]. Non-analytic effects from 
vector mesons are automatically taken into ac- 
count, because these vector qq bound states ap- 
pear as poles in the quark-photon vertex solu- 
tion [10]. 

We employ the model described in Sec. |], 
without any parameter adjustment. Since we work 
in the SU(2) isospin limit, u and d quarks are 
identical apart from their charge, and we have 
for the pion electromagnetic form factor simply 
F 7r {Q 2 ) = F uau (Q 2 ). Thus there are only three 
independent form factors, F U u U (Q 2 ), F uSu (Q 2 ), 
and F uSS (Q 2 ). With this model we have demon- 
strated explicitly that the impulse approximation 
is indeed independent of the unphysical parame- 
ter r\ and satisfies current conservation, provided 
that complete representations of the Dirac matrix 
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Figure 4: The pion charge form factor F n (Q 2 ) as ob- 
tained from different treatments of the quark-photon 
vertex. The inset shows the Q 2 region relevant for 
the charge radius. The data correspond to \F K \, taken 
from Refs. [26] (circles), [27] (squares), and [28] (dots). 



structure and the k • P dependence are used [11]. 
However, with the leading terms in the Chebyshev 
expansion of the BSAs and the quark-photon ver- 



tex given in Eq. (2.5) only, current conservation 
is violated, and the results do depend on the mo- 
mentum routing in the loop integral. Including 
terms up to order (k ■ P) 3 in this expansion re- 
stores current conservation and the independence 
of the momentum routing to within the numeri- 
cal accuracy which we estimate to be around 1%. 
Higher order terms do not change the results more 
than 1%, at least not at the values of Q 2 we have 
considered. Also a truncation of the BSAs or the 
quark-photon vertex in terms of their Dirac am- 
plitudes (E, F, G, and H for the pseudoscalar 
mesons, and eight transverse amplitudes for the 
quark-photon vertex) generally leads to a viola- 
tion of current conservation and/or independence 
of the choice of integration variable. 

In Fig. ^ we show our results for the pion form 
factor [10] around Q 2 = 0. With dressed prop- 
agators, the use of a bare quark-photon vertex 
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in Eq. (4.4) clearly violates electromagnetic cur- 
rent conservation and leads to F 7T (Q) ^ 1. Use of 
the Ball-Chiu Ansatz [29] for the dressed quark- 



Figure 5: Our results for the pion charge form fac- 
tor [11] compared to new data from JLab [30] including 
a re-analysis of older data. 



photon vertex conserves the electromagnetic cur- 
rent and ensures F(Q 2 = 0) = 1. However, the be- 
havior of the form factor away from Q 2 = is not 
constrained by current conservation, and in the 
present model, use of the Ball-Chiu Ansatz leads 
to a value for r 2 which is about 50% too small [10]. 
With the quark-photon vertex produced by solu- 
tion of the ladder BSE, and with quark propaga- 
tors from the rainbow DSE, all constraints from 
current conservation are satisfied and the calcu- 
lated value of r 2 is within 5% of the experimen- 
tal value [10]. Since the impulse approximation 
omits meson loop mechanisms such as pion rescat- 
tering, previously estimated to add up to 15% to 
the charge radius [25], the precision of the agree- 
ment shown in Fig. ^ should be considered prelim- 
inary until these aspects are investigated within 
the present model. 

In Fig. ||we show our results for Q 2 F 7r (Q 2 ) on 
a larger <5 2 -domain [11]. We obtain an excellent 
description of the spacelike Q 2 data, including the 
latest JLab results, without any adjustment of the 
model parameters. Our result on this Q 2 domain 
can quite well be described by a simple monopole, 
with a mass slightly lower than m p . 

The obtained charged kaon form factor is also 
in good agreement with the available data, see 
Fig. ^. The behavior of Fk on a larger space- 
like Q 2 -domain anticipated for future JLab [33] 
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0.08 0.12 
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Figure 6: The calculated K + form factor compared 
to the data from Refs. [31] (open diamonds) and [32] 
(solid squares). Within numerical errors, F K +(0) = 1. 



radii 


experiment 


calculated 


r 2 


0.44 ± 0.01 fm 2 


0.45 fm 2 


r 2 

' K+ 


0.34 ±0.05 fm 2 


0.38 fm 2 


r 2 


-0.054 ± 0.026 fm 2 


-0.086 fm 2 



Table 2: Our results for the charge radii, compared 
with the experimental values given in Refs. [28,32,34]. 

data, is discussed elsewhere [11]. The obtained 
charge radii for the pion and both the neutral and 
charged kaons are presented in Table ^, and agree 
well with the experimental data. These charge 
radii are somewhat larger than those obtained in 
a previous study [35] that was framed in terms 
of semi-phenomenological representations for the 
confined quark propagators and the BSAs within 
the impulse approximation. The main difference 
with that work is that here we use numerical solu- 
tions of truncated DSEs for all the elements needed 



in Eq. (4T), and that all our parameters were fixed 
previously. 

4.5 Timelike behavior 

Experimentally, F- K {Q 2 ) shows a resonance peak 



at Q 2 



-m 2 (in Euclidean metric, Q < cor- 



responds to the timelikc region). In general, for 
timclikc photon momenta Q 2 in the vicinity of a 
vector meson mass-shell, the pseudoscalar meson 
charge form factor Fp(Q 2 ) will exhibit a resonance 



peak associated with the propagation of interme- 
diate state vector mesons followed by the decay 
V -> PP. That is [10] 



F P (Q 2 ) 



g V pp m v 



9v (Q 2 + m 2 / — imy IV) 



(4.13) 



Here m v j gy is the V-j coupling strength fixed by 
the V — > e + e~ decay, gypp is the coupling con- 
stant for the V — > PP decay, and Ty is the V me- 
son width, which is dominated by the latter pro- 
cess. 

This vector meson mechanism in the charge 
form factors of the pion and kaon is generated by 
the pole structure in the quark-photon vertex ex- 



hibited in Eq. (4.7); the relationship of the elec- 



troweak decay constant fy in Eq. (4.7) with gy is 
f p m p = \/2m 2 /g p for the p°, and f^m^ = 3 rn^/g^ 
for the 4> [12]. In the present study, the pole is at a 
real timelike value of Q 2 since the ladder truncat- 
ion of the BSE kernel does not generate a width for 
any meson bound state. For example, one would 
have to supplement the ladder BSE kernel with 
the intermediate state 7T7t production mechanism 
to produce a width for the corresponding p-meson 
resonance in the vertex beyond the threshold for 
pion production, Q 2 < —4 m 2 in the timelike re- 
gion. 

Our results shown in Fig. [?] indeed indicate 
poles at Q 2 — —0.55 GeV 2 in the two form fac- 
tors involving u"fu coupling and a pole at Q 2 = 
— 1.1 GeV 2 in the form factor involving 57s cou- 




Figure 7: Timelike behavior of the individual form 
factors that make up the pion and kaon charge form 
factors. 
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experiment 


theory 




6.02 ± .05 


5.4 


9<t>KK 


4.64 ± .14 


4.3 



Table 3: Coupling constants for vector meson strong 
decays extracted from timelike charge form factors. 

pling. These are the p and (j) respectively and 
the coupling constants g p7r7r and g^KK can be ex- 



tracted from a fit of Eq. (4.13) to the calculated 
timelike charge form factors. The results are given 
in Table 0, and agree reasonably well with exper- 
iment. The theoretical or experimental width is 
a significantly greater fraction of the p mass than 
it is of the <f> mass. Thus the omission of meson 
loops in the BSE kernel can be expected to lead 
to greater error in the p BS amplitude than in 
the (f> BS amplitude. For this reason we speculate 
that this is also the likely explanation of why the 
present impulse approximation for gp^ is further 
from experiment than the corresponding result for 
g<pKK ■ 

4.6 The 7*7r/3 and 7*717 form factors 

The 7T° — > 77 decay rate is dictated by the axial 
anomaly as prescribed by electromagnetic gauge 
invariance and chiral symmetry. The resulting in- 
variant amplitude, for photon helicities Ai and A2, 
is 

2 i a r 



M 



Ai A2 



1.9/r 



x ^ a pe^ei x ^k la k 20 , (4.14) 

and ki are the polarization vector and 
momentum of the i th photon, and a cm is the 
electromagnetic fine structure constant e 2 /47r. As 



where e^^ 



indicated in Eq. (4.14), this result is convention- 
ally expressed in terms of / w = f n /y/2, the pion 
decay constant in the convention where its value 
is 92 MeV, rather than the convention /„■ = 131 
MeV used throughout this present work. The cou- 
pling constant g^-y-y appearing in Eq. (4.14) has 
the chiral limit value 3^J 77 = 1/2 and the resulting 
invariant amplitude for 7r°77 coupling saturates 
the axial anomaly through the massless pion pole 
in the divergence of the axial current as coupled 



to the electromagnetic field. The resulting decay 
width 



(q° ) 2 a 2 



167T3/2 



(4.15) 



is within 2% of the experimental width of 7.8 eV. 
The corrections due to finite pion mass are small. 

It is known that these features of the axial 
anomaly are preserved by the ladder-rainbow trun- 
cation of the DSEs combined with an impulse ap- 
proximation for the 7r°77 vertex because the rel- 
evant manifestations of electromagnetic gauge in- 
variance and chiral symmetry are present [24,36]. 

It is convenient to separate out the photon 
polarization vectors so that we deal with the 
vertex KJZ 1 defined by M XlX * = e^KJ^e^. 
Corresponding to the flavor decomposition 
7T° = (uu — dd) I \/2, the impulse approximation 
or triangle diagram for the vertex has the flavor 
decomposition 

A2J 7 (fc 1 ;fc 2 )= (Q 2 u A^-Q 2 d A d ^)/V2, (4.16) 

where Q u /d are the quark charges. We use isospin 
symmetry where A™ v = A^„. For each of these ver- 
tices there are two permutations of the final state 

photons to be included. The resulting expression 
; 2 



V2N C 



tr s [S(q)T*(q,k)S(k) 



x iTl(k,p)S(p)irZ(p,q)] . (4.17) 

Here the momenta p, k are determined in terms 
of the integration momentum q and the outgoing 
external momenta k\ , k 2 by momentum conserva- 
tion. 

To conveniently describe the 7*^7 transition 
and define the associated form factor, we replace 
k\ by the virtual photon momentum Q and express 
ki as &2 = P — Q/2 where the incoming pion mo- 
mentum is P + Q/2. Thus the new independent 
external momenta are P and Q and this facilitates 



The prefactor in Eq. (4.17) is different than that in 



Refs. [37,38] because a) the BSE normalization conventions 
are different and b) the symmetry factor of 2 from the iden- 
tity of the decay photons is here included in the vertex. 
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2.0 
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Figure 8: The 7* ^7 form factor [38], with data from 
CLEO and CELLO [40]. 

comparisons with our descriptions of meson elastic 
vertices and form factors. With one virtual pho- 
ton, the general form of the vertex is 



(4.18) 



where the defined form factor satisfies F(Q) = 1. 

The present model produces the coupling con- 
stant shown in Table |] at the physical pion mass; 
the chiral limit value is essentially the same within 
numerical accuracy. Our result [38,39] for the 
transition form factor is displayed in Fig. ||. The 
corresponding transition radius, shown in Table ^, 
agrees well with the experimental value reported 
by the CELLO collaboration [40]. For the mod- 
erate range of Q 2 displayed in Fig. [|, both the 
data and our DSE model are close to the monopole 
shape fitted to the Brodsky-Lepage asymptotic 
limit [41] Q 2 F(Q 2 ) ~ 4tt 2 /2 obtained from pQCD 
in the factorized approximation 3 . 

For the radiative decay p + — > 7r + 7, the photon 
can radiate from the u-quark or the ci-quark giving 

A^ = Q u A$« + Q d - , (4. 19) 

where A^' 9 is the vertex having the indicated 
quark flavor labelling and containing no charge 

3 This has been converted to the convention f-^ = 131 
MeV used in the present work. 



or flavor weights associated with external bosons. 
With both the p° and 7T° given by (uu- dd)/V2, 
the radiative decay p° — > 7r°7 can be expressed as 



A p *i = Q U A 



uu,u 



Q d Af/. 



(4.20) 



where the radiative contribution from both quark 
and antiquark of the same flavor have been com- 
bined. In the isospin symmetric limit, we have 
A uu,u = A dd,d = A ud,u = _A«^, and thus p ± 

and p° have identical 7r7 radiative decays at this 
level. In impulse approximation the physical ver- 
tex is [37] 

KV(P,Q) = ^ | A tr s [S(q)T*( q , q+ )S( q+ ) 

xiTl(q + ,q_)S(q-)TP(q^q)] , (4.21) 

where Q is the photon momentum and the p mo- 
mentum is represented by P — Q/2. This is com- 



pletely analogous to Eq. ( |4. L7[ ) for the 7*^7 vertex 
if the on-shell quark-photon vertex is replaced by 
the p BSA. In fact, if that photon momentum in 



Eq. (4.17) is continued into the timelike region, an 
analysis of the p° pole contribution in the manner 



discussed in Sec. 4.5 will identify Eq. (4.21) as the 
7*7r j o vertex. 

The p — > 7T7 coupling constant and the 7*7rp 
form factor are identified from the pit"/ vertex ac- 
cording to 

A^ViP, Q) = PcQpF^^iQ 2 ) . 

m p 

(4.22) 

At the photon point we have F(0) = 1 and g p7rj is 
the conventional coupling constant associated with 
the radiative decay width 



C^cm 9p-rr~y 

24 



TO p I 1 



(4.23) 



The radiative decay ui — * 7r°7 is given by 
A^{K, k 7 ) = Q u Ap u - Q d Afj d , (4.24) 



where, compared to Eq. (4.20), the change in 
phase of the second term here comes from the 
phase of the dd component of the ui. With 
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experiment 


theory 




0.501 ±.018 


0.50 


r 2 
7777 


0.42 ± .04 fm 2 


0.39 fm 2 


T-tt — >77 


7.8 ± 0.6 eV 


7.8 eV 




0.74 ± .05 GeV" 1 


0.69 GeV- 1 


1 p+ — >7T+7 


68 ± 7 keV 


53 keV 




102 ± 25 keV 


53 keV 


Quj7r~/ / ifbu) 


2.32 ± .08 GeV- 1 


2.07 GeV- 1 


^ U — >7T°7 


717 ± 43 keV 


479 keV 



Table 4: Coupling constant, interaction radius and 
width for the 7r° — > 77 decay along with coupling con- 
stants and widths for the indicated vector meson ra- 
diative decays. 



isospin symmetry, we have A^J 7 = A 
hence g^/m^ = 3g p7r7 /m p . 



P" and 



As Eq. (4.22) shows, it is g/my that is the 
natural outcome of the theory for the vector ra- 
diative decays, and it is this combination that we 
report in Table |[ We also give the corresponding 
decay widths. Except for the difference between 
the decay of the neutral and charged states of the 
p, which is beyond the reach of the isospin sym- 
metric impulse approximation, the agreement be- 
tween theory and experiment is within 10%. This 
is consistent with the other ir and p observables 
obtained from the same model. 

The 7*7T/9 form factor plays a role in the in- 
terpretation of electron scattering data from light 
nuclei, because the isoscalar meson-exchange cur- 
rent contributes significantly to these processes. 
In particular our understanding of the deuteron 
EM structure functions for Q 2 ss 2 — 6 GeV 2 re- 
quires knowledge of this form factor [42]. An 
initial exploratory study [43] of the j*np vertex 
within the present framework, but employing phe- 
nomenology to a much greater extent, produced 
a very soft result for the form factor and this is 
shown in Fig. ^ as the Mitchell-Tandy curve. It 
was found [42] that the resulting meson exchange 
current contribution provided a very good descrip- 
tion of the elastic deuteron electromagnetic form 
factors A(Q 2 ) and B(Q 2 ) in the range 2-6 GeV 2 
where such effects are important. 

The model under consideration here, based on 



> 

a 0.3 



0.2 



0.1 



current DSE calc. 
Mitchell-Tandy 1995 
VMD p monopole 




12 3 4 

Q 2 [GeV 2 ] 

Figure 9: The 7*7rp form factor. 

the DSEs of QCD, has no free parameter other 
than the two set by f v and (qq) as described ear- 
lier; the amount of phenomenology is significantly 
less than the earlier Mitchell-Tandy result [43]. 
The present work produces a form factor that is 
much softer than what is inferred from vector me- 
son dominance (VMD) but obviously harder than 
the Mitchell-Tandy result [43], as can be seen in 
Fig. ^. We expect the present impulse approxima- 
tion results for the form factor to be as reliable as 
the pion and kaon charge form factors from this 
model [11]. The influence of these recent results 
upon the analysis of electron scattering from the 
deuteron remains to be determined. In this regard, 
some estimation of the likely uncertainty from the 
virtuality of the ir and p within the assumed meson 
exchange current mechanism needs to be made. 

5. Extension to finite temperature 

The DSE approach can also be used to study QCD 
at finite temperature. This is done most easily 
via the Matsubara formalism where the energy 
components of the 4-momenta become discrete. 
Thus the quark propagator becomes a function of 
Pn = {p,u n ), where p is a three- vector, and u> n = 
(2n + 1)ttT is the fermion Matsubara frequency. 
The most general decomposition of the dressed 
renormalised quark propagator now becomes 

S^ 1 (Pn) = il -pA(p n ) + ij4UJ n C(p n ) + B(p n ) , 
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(5.1) 



which satisfies the DSE 



S 1 (p n ) = Z^ij-p + Z^ ij4UJ n + Z 4 m(n) 

+ S(P»), (5-2) 

with quark self-energy given by 



where we have employed the notation 



A oo , A ^ 



9./ 



T E 



(2tt) 3 



(5.4) 



Note that the argument of the gluon propagator, 
k n -i =Pn~Ql) corresponds to an even (bosonic) 
Matsubara frequency. For renormalization of the 
propagator we require 



S-\ P o 



17 -p + 174 lu + m(n) . (5.5) 



As before, we use the rainbow trunctation for the 
vertex, Y v = 7„, in combination with the ladder 
truncation for the BSE kernel K. 

5.1 Chiral symmetry restoration 

Chiral symmetry restoration has been studied at 
both zero and finite- (T, (j,) in a class of confining 
DSE models similar to the one discussed in Sees. || 
and [| [5]. Here, two models are considered in 
which the long-range part of the interaction is an 
integrable infrared singularity [44] motivated by 
T = studies [45] of the gluon DSE. Both models 
can be described by [13] 

9 2 D^(k m ) = P^V(k m ;m g ) + P^V(k m ;0), 



V(k rn ;m g ) = 2n 2 D^-5 0j 5 3 (k) 
+ V M (k* n + m 2 q ), 



(5.6) 



where P& = Pj t = 0, = % - hkj/k 2 , P£ v = 
Sfiv — k^k^/k^ — PT V , and m g is a Debye mass. 
The strength parameter D is fixed by a fit to m n 
and /„■ at T = 0. The models differ in the details 
of V M - 



It is anticipated that the restoration of chiral 
symmetry, which accompanies the formation of a 
quark-gluon plasma at nonzero temperature T, is 
a second-order phase transition in QCD with 2 
light flavors. Such transitions are characterised 
by two critical exponents: ([3,5), which describe 
the response of the chiral order parameters, X , 
to changes in T and in the current-quark mass, 
m. Denoting the critical temperature by T c , and 
introducing the reduced-temperature t := T/T c — 1 
and reduced mass h := m/T, then 



Xoc(-tf, t- 
X oc h 1/s , h 



0" , h = 0, 
> 0+ , 4 = 0. 



(5.7) 
(5.8) 



Calculating the critical exponents is an important 
goal because of the notion of universality, which 
states that their values depend only on the symme- 
tries and dimensions, but not on the microscopic 
details of the theory. 

One order parameter for dynamical chiral 
symmetry breaking is the quark condensate [3] 
(qq)^. There are other, equivalent order param- 
eters such as 

X ;=B(p = 0, wo), Xc ■= -Ty n V 5 ' 9 ) 

C(p = 0,LU O ) 

They should be equivalent, and the onset of that 
equivalence is a good way to determine the domain 
of reduced mass and reduced temperature in which 
the scaling behavior and critical exponents become 
manifest. In Ref. [13] it has been verified numeri- 
cally that in the chiral limit (m = 0) and for t ~ 0: 
/ w oc (qq) oc X(t, 0); i.e., that these quantities are 
all indeed bona fide order parameters. 

The success of the nonlinear a-model in de- 
scribing long-wavelength pion dynamics underlies 
a conjecture [46] that chiral symmetry restoration 
at finite T in 2-flavor QCD is in the same universal- 
ity class as the 3-dimensional, N — 4 Heisenberg 
magnet (0(4) model), with critical exponents [47] 
[3 H = 0.38 ± .01, 5 H = 4.82 ± .05. This is to 
be contrasted with mean-field critical exponents, 
/3 mf = 0.5, 5 mi = 3. 

5.2 Infrared dominant model 

A simple model [44] that is considered in re- 
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cent work [13,48-50] is the infrared dominant [ID] 
model with T>m(s) = in Eq. (5.6), with the mass- 



scale £> = 0.56 GeV fixed by fitting tt- and p- 
meson masses at T = 0. A current-quark mass 
of m = 12MeV yields m K = 140 Me V. There 
are two chiral limit solutions corresponding to 
two phases: the symmetric Wigner-Weyl solu- 
tion in which Bo(p n ) — 0, and the chirally broken 
Nambu-Goldstone solution in which 



2 D 



0, 



otherwise. 



From Eq. (p.7|), the critical exponent (3 for chiral 
restoration is simply (3 — 1/2 and T c — \/2D/2it 
~ 170 MeV. With a finite current-quark mass 
the quark DSE becomes a fourth-order algebraic 
equation and it is straightforward to establish [50] 
from Eq. ( |5.8| ) that 8 = 3. Clearly, this model has 
mean field critical exponents. The question is: is 
this result model-dependent or is it a more general 
property? 

5.3 1-loop perturbative QCD for the UV 

To provide a more realistic ultraviolet behavior, 
one can add a short-range term to the <5-function, 
corresponding to the exchange of a perturbative 
gluon at momenta p 2 ^> Aq CD . In Ref. [13], chi- 
ral symmetry restoration at finite temperature is 
studied using T>m{s) — G{s)/s in Eq. (5.6) where 



Q is given by Eq. (B.3) 4 . Again, the parameters 
are fixed at T = by requiring a good fit to m^iK 
and fn/ K [3]. 

The Debye mass is m 2 g = (16/5)w 2 T 2 and two 
values of the parameter oj are used which are phe- 
nomenologically equivalent at T — 0: u>x = 0.6m t 
and 102 = 1.2m t . To fix m n = 140 MeV, the former 
requires a renormalization point invariant current- 



4 Here a combination of a Gaussian function and <5- 
function is used to represent the strength of the interac tion 
at small momenta, whereas in Sec. Q we used Eq. ( |j.3| ) 
without an additional 5-function. Hence the results de- 
scribed here use D = 0.78 GeV 2 which is lower than the 
value D = 0.93 GeV 2 of the previous section. Qualitatively, 
the results are independent of the details of the effective 
interaction, as long as it is strong enough to break chiral 
symmetry at T = 0. 




110 120 130 140 150 160 

T (MeV) 



Figure 10: The analysis of 3 equivalent order param- 
eters for the critical exponent /3 within model u>2- 



quark mass of rh u = 6.6 MeV, while the latter re- 
quires rh u — 5.7 MeV. The critical temperature 
can be found from a calculation of the order pa- 
rameter in the chiral limit, rh — 0, as is shown in 
Fig. [n]. Once T c is known, one can analyse the 
behavior of the order parameter as a function of 
the current mass at T = T c , see Fig. [ll]. 

Although one can in principle extract the crit- 
ical exponents (3 and 8 from these plots, a much 
better method is to use the chiral and thermal 
susceptibilities [50]. The critical exponents ob- 
tained from the susceptibilities [13] are in good 
agreement with a direct calculation of the criti- 



cal exponents using Eqs. (5.7) and (5.8), and are 




Figure 11: The analysis of 2 equivalent order param- 
eters for the critical exponent 8 within model u>2- 
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ID 




U>2 


T c (MeV) 


169 


120 


152 


a 


i 

2 


0.50 


0.50 


5 


3 


3.00 


3.02 



Table 5: Critical temperature for chiral symmetry 
restoration and critical exponents characterising the 
second-order transition in the three illustrative mod- 
els [13]. 

given in Table [| Note that the scaling relations 
are only valid for very small current-quark masses: 
\og w {m/rh u ) < —5, where rh u corresponds to the 
current-quark mass that gives = 140 MeV in 
this model. It is only in this scaling region that 
the different order parameters become equivalent. 
These results are qualitatively, and for the criti- 
cal exponents quantitatively, independent of the 
parameters in this model. 



6. Spatial qq modes 

At T = the mass-shell condition for a meson as 
a qq bound state of the BSE is equivalent to the 
appearance of a pole in the qq scattering amp- 
litude as a function of P 2 . At T ^ in the 
Matsubara formalism, the 0(4) symmetry is bro- 
ken by the heat bath and one has P — > (P, fl m ) 
where fi m = 2mirT. Bound states and the poles 
they generate in propagators may be investigated 
through polarization tensors, correlators or Bethe- 
Salpeter eigenvalues. This pole structure is char- 
acterized by information at discrete points fl m on 
the imaginary energy axis and at a continuum of 
3-momcnta. Analytic continuation for construc- 
tion of real-time Green's functions (and related 
propagation properties) has been well-studied [51]. 
An unambiguous result is obtained by the require- 
ment that the continuation yield a function that 
is bounded at complex infinity and analytic off 
the real axis [51]. One may search for poles as a 
function of P 2 thus identifying the so-called spa- 
tial or screening masses for each Matsubara mode. 
These serve as one particular characterization of 
the propagator and the T > bound states. Here 



we discuss the mass dependence of the poles in the 
zeroth Matsubara mode, tt m = 0. 

6.1 Chiral modes near T c 

In the Matsubara formalism, the number of cou- 
pled equations in the the quark DSE, Eq. (ETj), 



and the meson BSE, Eq. (2.2), scales up with 
the number of fermion Matsubara frequencies in- 
cluded. At low temperature, the necessary large 
number of such modes combined with the loss of 
0(4) symmetry entails that a straightforward ex- 
tension of large scale T — calculations is pro- 
hibitive. Meson qq modes at T > have of- 
ten been studied within the Nambu— Jona-Lasinio 
model where the contact nature of the effective 
interaction allows decoupling of the Matsubara 
modes and analytic summation methods, see for 
example, Ref. [52] and references therein. How- 
ever the lack of quark confinement in that model 
leads to unphysical thresholds for qq dissociation. 
To learn about the fate of meson modes as influ- 
enced by deconfinement and chiral restoration, it is 
desirable to explore the T ^ properties of mod- 
els that, at T = 0, confine quarks and provide a 
realistic description of the light mesons. 

Such a study has recently been carried out [14] 
within the model of Sec. 5.3. This work was re- 
stricted to temperatures near and above the transi- 
tion, T > 100 MeV, where a maximum of 10 Mat- 
subara modes were found to be sufficient; only 
pseudoscalar and scalar modes were considered. 
In the chiral limit and at T = 0, there is dyman- 
ical chiral symmetry breaking and thus one finds 
massless pseudoscalar Goldstone bosons. The low- 
est scalar bound state in this model is an ideal- 
ized CT-meson with m a = 0.56 GeV. Such a low- 
mass scalar is typical of the rainbow-ladder trun- 
cation, although there is some model sensitivity. 
The rainbow-ladder truncation yields degenerate 
isoscalar and isovector bound states, and ideal fla- 
vor mixing in the 3-flavor case; improvements be- 
yond the ladder truncation are required in order to 
describe the observed scalar mesons. Furthermore, 
in the isoscalar-scalar channel, it will be necessary 
to include couplings to the dominant irir mode [53] 
because of the large a width. In the absence of 
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such corrections, the a properties discussed below 
are those of an idealized chiral partner of the tt. 

The poles in the Oth Matsubara modes of the 
scalar and pseudoscalar vertices, Ts(p n , k n ; P, 0) 
and rps(Pn, k n \ P, 0) where P = ft — k, evolve 
with T such that at the critical temperature, one 
has degenerate, massless scalar and pseudoscalar 
bound states [14]. Above T c , these bound states 
persist, becoming increasingly massive with in- 
creasing T. Furthermore, all but the leading Dirac 
amplitude vanish above T c and the surviving pseu- 
doscalar amplitude E v is pointwise identical to 
the surviving scalar one. These results indicate 
that the chiral partners are locally identical above 
T c , they do not just have the same mass. Ap- 
proaching T c from below, the a mass behaves as 
m a oc (1 — T/T C )P within numerical errors. Note 
that this is the same behavior as the chiral con- 
densate and other, equivalent order parameters 
for chiral symmetry restoration, such as the decay 
constant f % , as discussed in the previous section. 
For nonzero current-quark masses chiral symme- 
try restoration is exhibited as a crossover rather 
than a phase transition, see Fig. 12. The a mass 



exhibits a dip around the phase transition. The 
scalar and pseudoscalar meson masses become in- 
distinguishable at T ~ 1.2 T c . 

With the BSAs obtained from the homo- 
geneous BSE, one can now calculate the T- 




Figure 12: The T-dependence of and f w along 
with the mass of the scalar chiral partner of the pion 
(adapted from Ref. [14]). The model used is u>2, see 
Table @ and Sec El 




0.1 0.12 0.14 0.16 0.18 0.2 

T [GeV] 



Figure 13: T-dependence of the coupling T(0) and 
the 7T° — > 77 width (adapted from Ref. [14]). 



dependence of other meson properties [14]. The 
decay of a pion into two photons is calculated in a 
similar fashion as at T = in the previous section. 
Using the Ball-Chiu Ansatz [29] for the photon 
vertices, one can analytically show that in the chi- 
ral limit the pseudoscalar piece (E n ) of the pion 
BSA [36] saturates the Abelian anomalous contri- 
bution to the divergence of the axial- vector vertex: 



371-77 = \- The tensor structure of Eq. (4.21) sur- 
vives at nonzero temperature 5 



Aj4(fel,fc2 



= (fcixM,T(0), (6.1) 

7T 



where k\ and ki are the photon momenta (or 
rather, the spatial components of the zeroth Mat- 
subara modes of the photons), and with 7~(0) = 
g ni7 F(0)/f v . With a suitable extension of the 
Ball-Chiu Ansatz to nonzero temperature one can 
determine the T dependence of the coupling con- 
stant 3^77, or rather, of g vll j fi v which is the more 
interesting quantity. 

The T-dependence of the 7r°77 coupling and 
the associated width is depicted in Fig. [l3| Clear 
in the figure is that T(0) vanishes at T c ; in fact, 
it vanishes with the same critical exponents as 
the other order parameters: f3 = 0.5, indicating 
mean field behavior [14]. An accurate calcula- 
tion of the critical exponent is possible because the 
Goldberger-Treiman- like relation [9] f°E^{p n ) = 
Bo(p n ) is satisfied for all T. It is therefore not nec- 

5 In this section we use the convention f w = 92 MeV. 



16 



Confinement @ ESI, May-July 2000 



Pieter Maris and Peter C. Tandy 



essary to solve the BSE for the chiral pion bound 
state amplitude E^. Thus, in the chiral limit, 
the coupling to the dominant decay channel closes 
for both charged and neutral pions. These fea- 
tures were anticipated in Ref. [46]. Further, the 
calculated T(0) is monotonically decreasing with 
T, supporting the perturbative analysis to order 
(T 2 /f 2 ) in R ef - [54]. For rh ^ both the coupling: 
g w o 77 // 7r , and the width exhibit the crossover with 
a slight enhancement in the width as T — > T c due 
to the increase in m ff . 

The isoscalar-scalar-7T7r strong coupling van- 
ishes at T c in the chiral limit, which can be traced 
to .Bchirai — > 0. For in ^ 0, the coupling reflects 
the crossover. For the present case of the ideal- 
ized cr-meson, there are no physical consequences 
because the width for the decay a — > tttt vanishes 
just below T c where the a meson mass falls below 
2m n and the phase space factor vanishes. 

6.2 Separable model 

For studies at low finite temperature, where con- 
tinuity with T = results is to be verified, the re- 
quired number of Matsubara modes can be of order 
~ 10 3 . This is rather cumbersome with the model 
discussed above, but can be done if one utilizes a 
separable Ansatz [15,56] for the interaction that 
can implement the essential qualitative features 
of DSE/BSE models at finite temperature. Sep- 
arable representations have previously been found 
capable of an efficient modeling of the effective qq 
interaction in the infrared domain for T — meson 
observables [57,58]. 

In Ref. [15] a simple separable interaction is 
used, which confines quarks at T = 0. The few pa- 
rameters are fixed by tt and p/u> properties. The 
approach is a simplification of one developed ear- 
lier [58] that was found to be quite successful for 
the light meson spectrum at T = 0. The confin- 
ing mechanism is an infrared enhancement in the 
quark-quark interaction that is strong enough to 
remove the possibility of a mass shell pole in the 
quark propagator for real p 2 . In this implementa- 
tion, it is particularly transparent that sufficiently 
high temperature will necessarily restore a quark 
mass-shell pole and there will be a dcconfincmcnt 











T 






--- A(p=0) 






— B (p=0) 






B(p=0) 

— - C(p=0) 











T [GeV] 

Figure 14: T-dependence of quark self-energy am- 
plitudes at p — (adapted from Ref. [15]). 

transition. This model also implements low tem- 
perature dynamical chiral symmetry breaking and 
it preserves the Goldstone theorem in that the 
generated tt is massless in the chiral limit. The 
solutions of the BSE for the tt and p modes are 
particularly simple and are used to study the In- 
dependence of the meson masses and decay con- 
stants in the presence of both deconfinement and 
chiral restoration mechanisms. 

Here we discuss the results obtained [15] with 
the rank-2 separable interaction 

Df v (p~q) = <WA) Mp 2 )fo(q 2 ) 

+ D 1 h( P 2 )(p-q)h(q 2 )] , (6.2) 

and a Feynman-like gauge is chosen for phe- 
nomenological simplicity. The choice for the two 
strength parameters Dq,Di, and corresponding 
form factors fi(p 2 ) is constrained by consideration 
of the resulting solution of the DSE for the quark 
propagator in the rainbow approximation, and fit- 
ted to reproduce n and p properties at T = 0. 
In Ref. [15] a single Gaussian representation for 
the fi(p 2 ) was used to represent the major fea- 
tures. The finite T extension is implemented us- 
ing the Matsubara formalism as before, via the dis- 
cretization of the timelikc component of momenta: 

The T-dependence of the solutions for A, B 



and C at p 1 = is displayed in Fig. 14. The re- 
sults shown are for the lowest Matsubara mode 
(ujq = irT) which provides the leading behavior 
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as T is increased. The chiral restoration critical 
temperature T c — 121 MeV is identified from the 
vanishing of the chiral limit amplitude Bq(Q,ujq) 
as shown. Below T c , A and C are relatively con- 
stant, 0(4) symmetry is approximately manifest, 
and the main effect is an almost constant quark 
wave function renormalization via 1/C; the cen- 
tral feature is a rapidly decreasing mass function 6 . 
Above T c , there remains a significant temperature 
range where the self-interaction effects are strong, 
both A and C are considerably enhanced above 
their perturbative values, and the breaking of 0(4) 
symmetry is manifest. The present model thus 
captures the qualitative T-dependence observed 
for the dressed quark propagator in studies of the 
quark DSE [50,59]. 

The deconfinement transition is characterized 
by the appearance of a quark progagator pole for 
real values oip 2 . The chiral restoration point T c 
and the quark deconfinement point T4 are gen- 
erally expected be to identical or nearly so [60]. 
The separable model under consideration here pro- 
duces T d = 0.9 T c = 105 MeV. 

6.3 Vector mesons at nonzero T 

Due to the breaking of 0(4) invariance, the T = 
transverse vector meson splits for T > into 3- 
space longitudinal and transverse modes. For the 
spatial modes characterized by (P, fin = 0), the 
BSAs studied in the separable model are [15] 

T^ L \ Pn ,k n -P) = 6^j 4 fo(q 2 n )F p{L) (P 2 ) , (6.3) 

and 

rf T \ Pn ,k n ;P) = 

U-^^)MQ 2 n)F p( T)(P 2 ), (6.4) 

where q n = (p n + k n )/2 is the relative momen- 
tum. The T-dependence of the corresponding 
masses is displayed in Fig. [l5| These modes 
are almost degenerate and T-independent until 

6 The rank-1 limit of this model has A = C = 1 for all T 
and the behavior of B(Q,wo) as function of T is similar to 
that in Fig. F3 except T c = 146 MeV is obtained [15] . 




0.4 - 
0.2 - 



0.05 0.1 0.15 0.2 

T [GeV] 

Figure 15: T-dependence of the longitudinal p mass 
(dot-dashed line) and the transverse p mass (solid 
line). Also shown is 2M W indicating that the strong 
decay channel becomes inaccessible within 20-30 MeV 
beyond T c (adapted from Ref. [15]). 

about T c /2 where the breaking of 0(4) invari- 
ance becomes significant. The qualitative fea- 
tures M^(T) > Mj(T) and Mj(T) « constant 
for T < T c seen here in the context of a finite range 
interaction have previously been noted within the 
limiting case of the zero momentum range ID 
model [49]. This latter model was not studied for 
T > T c . In the separable model discussed here, 
the longitudinal mode becomes unstable to qq dis- 
sociation at T ~ 180 MeV; the transverse mode 
continues to be below the spatial qq threshold for 
the temperature range displayed. 

Next, one can study the couplings appropri- 
ate to p — * e + and p — * irir for the meson 
modes characterized by P = (P, 0) [15]. The var- 
ious spatial modes from fi m 7^ are characterized 
by masses much greater than that of the lowest 
mode considered here. It is anticipated that this 
lowest mode characterizes the qualitative behav- 
ior of the physical decay processes p° — > e+ e~ and 

p — * TTTT. 

The T — expression for the electromagnetic 
coupling constant g p has a straightforward exten- 
sion to T > for a transverse spatial p° mode [15]. 
Use of the f2 m = solution described above yields 

Mf{T) 
9p(T) 
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leads to the electromagnetic decay width 



T [GeV] 



Figure 16: T-dependence of the transverse p par- 
tial widths due to electromagnetic e + e~ decay (dashed 
line) and strong nn decay (solid line) (adapted from 
Ref. [15]). 



Nr. 



tr s [74 S(q n+ ) rf T \q n+ ,q n -;P) S(q n -)] . 

(6.5) 



The impulse approximation for the pun ver- 
tex [4,37], after extension to T > for spatial 
modes characterized by Q = (Q, 0) for the p and 
P = (P, 0) for the relative 7T7t momentum, takes 
the form 

A l/ (P,Q)=P u9pw7r (T) = 

-2N C / tr s [S(q) T n (q, q n+ ; -P+) S(q n+ ) 

x K {T) (<Zn+ ,q n -; Q) %»- ) ( Qn _ , q; P_j\ , 

(6.6) 

where P± — P ± Q/2. Use of Eq. ( |6~3| ) immedi- 
ately shows that the longitudinal p mode cannot 
couple to 7T7T. With summation over enough Mat- 
subara modes for convergence, both the electro- 
magnetic coupling g p (T), and the strong coupling 
gp-KTriT), reproduce the independently determined 
T = results of this model. There is very little In- 
dependence below about 0.8 T c where there is ap- 
proximate 0(4) symmetry as evident in the quark 
propagator behavior in Fig. Near T c , g p (T) 
rises and gp^fT) decreases. 

Combined with the relevant phase space, this 



r p o_ 



eAT) 



(6.7) 



while the corresponding strong decay width is 



r(T) 



g% n (T)Mj(T) 



4tt 



12 



1 - 



Mf{T) 



-i 3/2 



(6.8) 

The T-dependence estimated in this way for the 
decay widths is due to the response of the quark 
substructure to the heat bath, particularly the res- 
toration of chiral symmetry. The results are dis- 
played in Fig. The contrast between the behav- 
ior of the electromagnetic and strong widths near 
and just above T c should be a more robust finding 
than the details of the individual processes. The 
strong width decreases rapidly and vanishes just 
above T c while the electromagnetic width remains 
within 20% of the T = value. Part of the strong 
decrease of the intrinsic 7T7t width of the transverse 
p is due to the decrease in the coupling constant, 
however the dominant effect is the T-dependence 
of the phase space factor in Eq. (6.8). As displayed 
in Fig. m, 2M,r(T) rises faster with T than does 
Mj(T) until at T = 1.17 T c one has Mj = 2AU. 
Beyond this point, the strong decay p T — > irir is 
phase-space blocked. This suggests that the total 
p T width of 151 MeV at T = decreases by about 
50% near T — T c and drops sharply to the electro- 
magnetic value of about 6 keV by T = 1.17 T c . 

One would expect the mass of the pseu- 
doscalar K correlation to rise with T in a similar 
fashion to M^, while the masses of the vector 4> and 
K* modes should rise like the p. This suggests that 
the vector modes p, K* and 4> tend to be trapped 
with their relatively long electroweak lifetimes and 
with significantly increased masses for a domain of 
high temperatures above the transition. This sug- 
gests that within the gas of pions and other pseu- 
doscalars that dominate the hot hadronic product 
from heavy-ion collisions, the role of vector me- 
son correlations in producing the dilepton spectra 
could be significantly less than conventional ex- 
pectations. 
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This narrowing of the intrinsic decay width of 
the vector meson mode in the heat bath, and also 
the decrease of the decay width of the idealized o- 
meson discussed in Sec. |6.1[ , is a mechanism that is 
distinct from the collisional broadening effect [61] 
from the many-hadron environment. The results 
discussed here indicate that there is a non-trivial 
T-dependence to intrinsic coupling constants, such 
as <7 p7rw and (fo^r, and decay phase space. The 
intrinsic effect tends to significantly decrease the 
decay width; the hadronic medium effects have the 
opposite influence. An approach that incorporates 
both phenomena is clearly called for. 

Also note that only the spatial or screening qq 
masses have been investigated within this model. 
The temporal masses provide a different charac- 
terization of the correlations; lattice simulations 
indicate that spatial masses become much larger 
than the temporal masses above T c [55] , whereas a 
Nambu-Jona-Lasinio model study [52] found that 
they are significantly different only in the range 
150 MeV < T < 350 MeV. 

6.4 Behavior at high T 

Above the deconfinement and chiral transition 
temperature, it might be expected that meson 
modes should dissolve in favor of a gas of essen- 
tially massless quarks. However for a significant 
temperature range above T c , the spatial ir and p 
modes studied here continue to be stable against 
qq dissociation and do not dissolve into a free 
quark gas [15]. The results for M„(T) and Mj(T) 



are displayed in Fig. 17, along with the quark dy- 
namical mass function at p — 0, using the simpler 
rank-1 separable interaction. The masses of both 
spatial meson modes approach the asymptotic be- 
havior 2ttT from below. This asymptotic behavior 
has been observed in lattice simulations [60,62] as 
well. 

The emergence of the 2irT behavior for masses 
of qq correlations at large T can easily be under- 
stood in the rank-1 model, as is shown in Ref. [15]. 
There it is found that for spatial meson modes 
M (T) ~ 2ttT- A(T) where A is a positive mass 
defect that will typically decrease with T. Thus 
the spatial meson mass or screening mass will ap- 
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Figure 17: High T behavior of the spatial masses of 
the 7r and transverse p modes from the rank-1 separa- 
ble model (adapted from Ref. [15]). 

proach the thermal mass of a pair of massless 
fermions from below. 

The same qualitative behavior at large T can 
be seen in a semi-analytic way by considering the 



ID model of Sec. 5.2. For example, it is found [15] 



that the leading large T behavior of Mj(T) is 



Mj\T) - (27TT) 2 (1- JL 



(6.9) 



Due to the long range nature of this model, the 
asymptotic spatial mass defect of the transverse p 
mode actually becomes constant A p (T) — > r). 

the 7r and transverse p spatial masses 



In Fig. 
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at T > T c from both the rank-1 separable model 
and the ID model are presented in comparison 
with lattice QCD simulations of spatial screening 
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Figure 18: Spatial masses from the rank-1 separable 
model (sep) and the ID model from Ref. [15] compared 
to lattice QCD simulations (latt) taken from Ref. [62] . 
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masses [62]. The solid horizontal line marks 2tt 
while the lower horizontal dot-dashed line repre- 
sents the lattice free limit corrected [62] for the 
lattice time extent N t . From Fig. [l^ it is evident 
that for the ir at T > T Cl the spatial or screening 
mass defect A — 2nT — M(T) is decreasing more 
rapidly above T c in the rank-1 separable model 
than is evident from the lattice simulations. This 
is consistent with the exponential behavior of the 
employed form factors. It is also consistent with 
the absence of quark vector self-energy amplitudes 
A{p n ) — 1 and C{p n ) — 1 through which interac- 
tions can be quite persistent in the asymptotic re- 
gion. This can be demonstrated within the chiral 
limit ID model where those amplitudes are strong 
and indeed have the power law fall-off. As seen 
in Fig. [l8|, the resulting mass defect for both it 
and p is in fact too strong when compared with 
the lattice QCD simulations. This persistent self- 
interaction well above T c , which slows the approach 
to free behavior such as Stefan-Boltzmann ther- 
modynamics [48] , may well be what is signaled by 
the lattice QCD data in Fig. [jj]. 

7. Summary 

Within a Poincare-covariant approach to modeling 
QCD based on the DSE framework, the descrip- 
tion of the ground state pseudoscalar and vector 
mesons has been reviewed at both zero temper- 
ature and finite temperature. The only approxi- 
mation made is a self-consistent truncation of the 
set of DSEs, which respects the relevant vector 
and axial-vector WTI. In particular we choose the 
truncation wherein the employed quark propaga- 
tors, the meson BSAs, and the quark-photon ver- 
tex are solutions of their DSEs in rainbow-ladder 
approximation. The two parameters are fixed pre- 
viously by fitting the chiral condensate, m^/K and 
fx. We include all relevant Dirac amplitudes for 
the BSAs and their full dependence upon the an- 
gular variable k ■ P. 

The obtained pion and kaon electromagnetic 
form factors are in good agreement with the avail- 
able data over the entire Q 2 range considered, and 
the calculated charge radii are within the bounds 



of the experimental values. The electromagnetic 
current is explicitly conserved in this approach, 
and there is no fine-tuning needed to obtain (0) 
= 1 = F K +(0) and F K o(0) = 0. We also demon- 
strate explicitly that our results are (within nu- 
merical accuracy) independent of the momentum 
partitioning of the BSAs. The timelike behavior 
of the charge form factors exhibits the vector me- 
son production resonances and this allows the as- 
sociated vector meson strong decay constants to 
be extracted. The confining property of the quark 
self-energy dynamically produced in this work pre- 
vents the appearance of spurious qq production ef- 
fects. We also discuss transition form factors for 
the processes j*p — > tt and 7*7r — > 7 as well as the 
coupling constants for the associated physical ra- 
diative decay of the p and the electromagnetic de- 
cay of 7T° as dominated by the axial anomaly. 

We discuss the finite temperature extension of 
the approach as implemented by the Matsubara 
method. Both the chiral restoration and the de- 
confinement transitions are identified and a study 
of the critical exponents for the chiral transition 
is reviewed. Several versions of the model are em- 
ployed to gain complementary viewpoints on the 
spatial qq correlative modes with (pseudo-) scalar 
and vector meson quantum numbers. The corre- 
sponding masses M„(T), Mj{T) and M%(T) are 
found to be almost T-independent below T c fol- 
lowed by a strong increase. This behavior is char- 
acteristic of lattice QCD simulations [60,62] and 
DSE studies [14,15,49]. The scalar mass M CT (T), 
an idealized chiral partner of the pion, decreases 
with T until it becomes equal to the pion mass 
at T c ; above the critical temperature for chiral 
symmetry restoration the scalar and pseudoscalar 
modes are degenerate. Once the scalar mass drops 
below twice the pion mass, the decay width as- 
sociated with the decay a — ► rnr drops to zero, 
because of phase-space blocking [14,63]. Also the 
strong decay p — ► tttt decreases with temperature, 
and drops to zero at about 25 MeV above T c due 
a phase-space blocking effect. 

One should keep in mind that the results pre- 
sented here follow from the response of the quark 
gluon content of the mesons to the heat bath. A 
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different phenomenon is the coupling of the meson 
modes to the many-hadron environment which in- 
troduces collisional broadening [61]. In order to 
study this effect, one has to go beyond rainbow- 
ladder truncation, and incorporate meson loop cor- 
rections. Also at zero temperature such correc- 
tions need to be studied in more detail: pion and 
kaon loops give rise to a nonzero width for the 
vector mesons. This would also allow for more re- 
alistic calculations of electromagnetic form factors 
in the timelike region, near and beyond the vector 
meson resonances. However, it is unclear how to 
add such corrections while preserving all relevant 
Ward identities. 
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